Prom a few to many electrons in quantum dots under strong magnetic fields: 
Properties of rotating electron molecules with multiple rings 
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Using the method of breaking of circular symmetry and of subsequent symmetry restoration via 
projection techiques, we present calculations for the ground-state energies and excitation spectra of 
iV-electron parabolic quantum dots in strong magnetic fields in the medium-size range 10 < N < 30. 
The physical picture suggested by our calculations is that of finite rotating electron molecules 
(REMs) comprising multiple rings, with the rings rotating independently of each other. An analytic 
expression for the energetics of such non-rigid multi-ring REMs is derived; it is applicable to arbitrary 
sizes given the corresponding equilibrium configuration of classical point charges. We show that the 
rotating electron molecules have a non-rigid (non-classical) rotational inertia exhibiting simultaneous 
crystalline correlations and liquid-like (non-rigidity) characteristics. This mixed phase appears in 
high magnetic fields and contrasts with the picture of a classical rigid Wigner crystal in the lowest 
Landau level. 

PACS numbers: 73.21. La; 71.45.Gm; 71.45.Lr 



I. INTRODUCTION 

A. Computational motivation 

Due to the growing interest in solid-state nanostruc- 
tures, driven by basic research and potential technologi- 
cal considerations, two-dimensional iV-electron semicon- 
ductor quantum dots (QDs) in field-free conditions and 
under applied magnetic fields (B) have been extensively 
studied in the last few years, both experimentally 1,2 ' 3 
and theoreticallyi^SiLSiS Experimentally, the case of 
parabolic QDs with a small number of electrons (N < 30) 
has attracted particular attention, as a result of pre- 
cise control of the number of electrons in the dot that 
has been demonstrated in several experimental investi- 
gations. 

Naturally, QDs with a small number of electrons are 
also most attractive for theoretical investigations, since 
their ground-state properties and excitation spectra can 
be analyzed^SiLSiiSiiiii^ through exact-diagonalization 
(EXD) solutions of the many-body Schrodinger equa- 
tion. In particular, in combination with certain approx- 
imate methods, which are less demanding computation- 
ally while providing highly accurate results and a trans- 
parent physical picture (e.g., the method of successive 
hierarchical approximations*^ see below), EXD calcula- 
tions confirmed the spontaneous formation of finite ro- 
tating electron molecules (REMs) and the description 
of the excited states with magic angular momenta as 
yrast rotational bands of these REMs^ (sometime the 
REMs are referred to as "rotating Wigner molecules," 
RWMs). However, the number of Slater determinants 
in the EXD wave-function expansion increases exponen- 
tially as a function of N, and as a result EXD calculations 
to date have been restricted to rather low values of N, 
typically with N < 10; this has prohibited investigation 
of REMs with multiple rings. A similar problem appears 
also with other wave functions that are expressed as a dis- 



crete sum over Slater determinants, such as the analytic 
REM wave functions [7[a,b)], or the variational Monte 
Carlo approach of Ref. UM 

Most EXD calculations (see, e.g., Refs. H l^b). ITflTlL 
Il2|) have been carried out in the regime of very strong 
magnetic field (i.e., B — ► oo), such that the Hilbert space 
can be restricted to the lowest Landau level (LLL): in 
this regime, the confinement docs not have any influ- 
ence on the composition of the microscopic many-body 
wave function (see section II. B). EXD calculations as a 
function of B that include explicitly the full effect of the 
confinement i.e., mixing with higher Landau lev- 
els are more involved, and thus they are scarce and are 
usually restricted to very small sizes with N < 4. An 
exception is presented by the method of hyperspherical 
harmonics?^ which, however, may not be reliable for all 
the sizes up to N ~ 10 (see below). 

Systematic EXD calculations beyond the numerical 
barrier of N ~ 10 electrons are not expected to become 
feasible in the near future. In this paper, we show that 
a microscopic numerical method, which was introduced 
by us recently and is based on successive hierarchical ap- 
proximations (with increasing accuracy after each step) is 
able to go beyond this barrier. This approach (referred 
to, for brevity, as the "two-step method") can provide 
high-quality calculations describing properties of QDs as 
a function of B in the whole size range 2 < N < 30, 
with (or without) consideration of the effect of the con- 
finement on the mixing with higher Landau levels. In 
this paper, we will consider the case of fully polarized 
electrons, which in typical GaAs experimental devices is 
appropriate for strong B such that the ground-state an- 
gular momentum L > Lq = N(N — l)/2 (see section II. A 
and footnotes therein). 

The minimum value Lq specifies the so-called maxi- 
mum density droplet (MDD); its name results from the 
fact that it was originally defined 15 in the LLL where 
it is a single Slater determinant built out of orbitals 
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with contiguous single-particle angular momenta 0,1,2, 
N—l. We will show, however, (see Section IV. B) that 
mixing with higher Landau levels is non-negligible for 
MDD ground states that are feasible in currently avail- 
able experimental quantum dots; in this case the electron 
density of the MDD is not constant, but it exhibits os- 
cillations. 



B. Nonclassical (non-rigid) rotational inertia 

The existence of an exotic supersolid crystalline 
phase with combined solid and superfluid characteris- 
tics has been long conjectured 1 ^ 7 ^ for solid 4 He un- 
der appropriate conditions. The recent experimental 
discoverjiiS that solid 4 He exhibits a nonclassical (non- 
rigid) rotational inertia (NCRI 18 ) has revived an intense 
interes*2fliSiiS2iSi24 in the existence and properties of the 
supersolid phase in this system, as well as in the possible 
emergence of exotic phases in other systems. 

As we show here, certain aspects of supersolid behavior 
(e.g., the simultaneous occurrence of crystalline correla- 
tions and non-rigidity under rotations) may be found for 
electrons in quantum dots. As aforementioned, under a 
high magnetic field, the electrons confined in a QD lo- 
calize at the vertices of concentric polygonal rings and 
form a rotating electron molecule^ We show that the 
corresponding rotational inertia strongly deviates from 
the rigid classical value, a fact that endows the REM 
with supersolid- like characteristics (in the sense of the ap- 
pearance of a non-classical rotational inertia, but without 
implying the presence of a superfluid component). Fur- 
thermore, the REM at high B can be naturally viewed as 
the precursor of a quantum crystal that develops in the 
lowest Landau level (LLL) in the thermodynamic limit. 
Due to the lack of rigidity, the LLL quantum crystal ex- 
hibits a "liquid" -like behavior. These conclusions were 
enabled by the development of an analytic expression for 
the excitation energies of the REM that permits calcu- 
lations for an arbitrary number of electrons, given the 
classical polygonal-ring structure in the QD. 25 

The paper is organized as follows. Section II is devoted 
to a description of computational methods for the proper- 
ties of electrons in QDs under high magnetic fields, with 
explicit consideration of effects due to the external con- 
finement. In section III, we compare results from various 
computationals methods with those obtained via exact 
diagonalization. Illustrative examples of the formation of 
crystalline rotating electron molecules with ground-state 
multiple concentric polygonal ring structures, and their 
isomers, are given in section IV for QDs with A = 6, 
9, 11, 17. The yrast band of rotational excitations (at a 
given B) is analyzed in section V along with the deriva- 
tion of an analytic formula that provides for stronger 
fields (and/or higher angular momenta) accurate predic- 
tions of the energies of REMs with arbitrary numbers of 
electrons. In section VI, we discuss the non-rigid (liquid- 
like) characteristics of electrons in quantum dots under 



high magnetic fields as portrayed by their non-classical 
rotational inertia. We summarize our findings in section 
VII. For an earlier shorter version of this paper, see Ref. 

M 

II. DESCRIPTION OF COMPUTATIONAL 
METHODS THAT CONSIDER THE EXTERNAL 
CONFINEMENT 

A. The REM microscopic method 

In our method of successive hierarchical approxima- 
tions, we begin with a static electron molecule (SEM), de- 
scribed by an unrestricted Hartree-Fock (UHF) determi- 
nant that violates the circular symmetry. 9 Subsequently, 
the rotation of the electron molecule is described by a 
post-Hartree-Fock step of restoration of the broken cir- 
cular symmetry via projection techniques. 7 9 Since we fo- 
cus here on the case of strong B, we can approximate the 
UHF orbitals (first step of our procedure) by (parameter 
free) displaced Gaussian functions; that is, for an electron 
localized at R 3 (Zj), we use the orbital 

1 / \z-ZA 2 , A f . 

u ^ z ' j > = "T^a exp I — 2A 1 M z i Z o; B )J j (!) 

with A = y/h/m*£l; f2 = where uj c — 

eB/(m*c) is the cyclotron frequency and loq specifies 
the external parabolic confinement. We have used com- 
plex numbers to represent the position variables, so that 
z = x + iy, Zj = Xj + iYj. The phase guaran- 
tees gauge invariance in the presence of a perpendicular 
magnetic field and is given in the symmetric gauge by 
<p(z, ZfB) = (xYj - yXj)/2l 2 B , with l B = y/he/eB. 

For an extended 2D system, the Zj's form a triangu- 
lar latticed For finite A, however, the Zj's coincid o 7 i 9 i 28 
with the equilibrium positions [forming r concentric reg- 
ular polygons denoted as (m, ri2, n r )] of A = J2q=i n q 
classical point charges inside an external parabolic 
confinement ~ In this notation, n\ corresponds to the 
innermost ring with ri\ > 0. For the case of a single 
polygonal ring, the notation (0, A) is often used; then it 
is to be understood that ni = N. 

The wave function of the static electron molecule 
(SEM) is a single Slater determinant |* SEM [z]) made 
out of the single-electron wave functions u{zi,Zi), i = 
1,...,N. Correlated many-body states with good total 
angular momenta L can be extracted 7 (second step) from 
the UHF determinant using projection operators. The 
projected rotating electron molecule state is given by 

I*?™) = / ... / d 7l ...d lr 

JO JO 

x|* SEM ( 7l ,..., 7r ))exp ^£ 7? ^j.(2) 

Here L = Y^ q =\ L q and |* SEM [7]) is the original Slater 
determinant with all the single- electron wave functions 
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of the qth ring rotated (collectively, i.e., coherently) by 
the same azimuthal angle 7(? . Note that Eq. @ can 
be written as a product of projection operators acting 
on the original Slater determinant [i.e., on |\l/ SEM (7i = 
0, ...,7 r = 0))]. Setting A = IbV% restricts the single- 
electron wave function in Eq. Q to be entirely in the 
lowest Landau levefi (see Appendix A). The continuous- 
configuration-interaction form of the projected wave 
functions [i.e., the linear superposition of determimants 
in Eq. (2)] implies a highly entangled state. We require 
here that B is sufficiently strong so that all the electrons 
are spin-polarized^ and that the ground-state angular 
momentum L > Lo = N(N — l)/2 (or equivalcntly that 
the fractional filling factor v = L/Lo < 1). 

Due to the point-group symmetries of each polygonal 
ring of electrons in the SEM wave function, the total 
angular momenta L of the rotating crystalline electron 
molecule are restricted to the so-called magic angular 
momenta, i.e., 

r 

L rn = L + ^ k q n qi ( 3 ) 

9=1 

where the k q 's are non- negative integers^ (when m = 1, 
ki = 0). 

The partial angular momenta associated with the qth 
ring, L q [see Eq. J2J], are given by 

L q = L 0i q + k q n q , (4) 

where L , 9 = Z)S™+i( i_ l ) witn *9 = Z)«=i n s (*i = °)j 
and L = YJq=i L o,q- 

The energy of the REM state [Eq. 10)] is givewLS by 

p27T I p27T 

E l EM =J o h{b]yWWd[l\ / J n([ 7 ])e lW - [L1 d[ 7 ], 

(5) 

with the hamiltonian and overlap matrix elements 
Kb}) = (* SEM ([0])|iJ|tf SEM ([ 7 ])) and n([ 7 ]) = 
(* SEM ([0])|* SEM ([ 7 ])), respectively, and [7] • [L] = 
J2q=ilqLq- The SEM energies are simply given by 
E S em = h([Q])/n([0]). 

The many-body Hamiltonian is 

N N N 2 

with 

being the single-particle part. The hamiltonian H de- 
scribes N electrons (interacting via a Coulomb repulsion) 
confined by a parabolic potential of frequency ujq and 
subjected to a perpendicular magnetic field B, whose 
vector potential is given in the symmetric gauge by 
A(r) = ^(—By, Bx, 0). m* is the effective electron mass, 



k is the dielectric constant of the semiconductor mate- 
rial, and rij = |rj — Tj\, For sufficiently high magnetic 
fields, the electrons are fully spin-polarized and the Zee- 
man term [not shown in Eq. 10] does not need to be 
considered^ Thus the calculations in this paper do not 
include the Zeeman contribution, which, however, can 
easily be added (for a fully polarized dot, the Zeeman 
contribution to the total energy is Ng* fisB /2, with g* 
being the effective Lande factor and ^jlb the Bohr mag- 
neton). 

The crystalline polygonal-ring arrangement 
(ni, ri2, n r ) of classical point charges is portrayed 
directly in the electron density of the broken-symmetry 
SEM, since the latter consists of humps centered at 
the localization sites Zj's (one hump for each electron). 
In contrast, the REM has good angular momentum 
and thus its electron density is circularly uniform. To 
probe the crystalline character of the REM, we use the 
conditional probability distribution (CPD) defined as 

P(r,r ) = ($|^(r 4 -r)5(r, - r )|$)/($|<f), (8) 

where $(ri, r 2 , rjy) denotes the many-body wave func- 
tion under consideration. F(r,r ) is proportional to the 
conditional probability of finding an electron at r, given 
that another electron is assumed at ro. This procedure 
subtracts the collective rotation of the electron molecule 
in the laboratory frame of referenece, and, as a result, 
the CPDs reveal the structure of the many body state in 
the intrinsic (rotating) reference frame. 

B. Exact diagonalization in the lowest Landau level 

We describe here a widely used approximationiii2i2i 
for calculating the ground state at a given B, which takes 
advantage of the simplifications at the B — > 00 limit, i.e., 
when the relevant Hilbert space can be restricted to the 
lowest Landau level [then huio << hu> c /2 (for B — > 00) 
and the confinement can be neglected at a first step]. 
Then, the many-body hamiltonian [see Eq. |JBJ] reduces 
to 

* N N 2 

( 9 ) 

2—1 J>i J 

Due to the form of the limiting Hamiltonian in Eq. JJjJ , 
one can overlook the zero-point-energy term and perform 
an exact diagonalization only for the Coulomb interac- 
tion part. The corresponding interaction energies can be 
written as 

E^(L)=S^(L)^, (10) 

kIb 

where E^t^LL ^ s dimensionless. The subscript "int" iden- 
tifies the e — e interaction as the source of this term. 
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In this approximation scheme, at finite B the external 
confinement hiog is taken into consideration only through 
the lifting of the single-particle degeneracy within the 
LLL, while disregarding higher Landau levels. As a re- 
sult, the effect of the confinement enters here only as 
follows: (I) in the interaction term [see Eq. QlOjl]. one 
scales the effective magnetic length, i.e., one replaces Ib 
by A/v2 (see section II. A for the definition of Is and A) 
without modifying the dimensionless part ^^lllj an< ^ 
(II) the contribution, E^=°{B,L) (referenced to NMl), 
of the single-particle hamiltonian 53 i=1 H sp (i) to the to- 
tal energy [see Eq. @J] is added to E^® LL (L) [corre- 
sponding to the second term on the right-hand side of 
Eq. ©]. E™=°{B,L) is the sum of Darwin-Fock single- 
particle energies e^f with zero nodes (n — 0; the cor- 
responding single-particle states become degenerate at 
B — > oo and form the lowest Landau level). Since 



DF 
C n,l. 



= (2n+l+ \l\)%n- ITiuj c /2, 



(11) 



the E™~°(B,L) is linear in the total angular momentum 



E?=°(B,L) = h(Sl-u; c /2)L. 



(12) 



Note that E™=°(B -> oo,L) -> 0. 

We denote the final expression of this approximation 

by -^to^LLL; it is given by 
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FIG. 1: (Color online) Two-step-method versus EXD calcula- 
tions: Ground-state energies for N = 4 electrons (referenced 
to 4ftSl) as a function of the magnetic field B. Thick dashed 
line (red): broken-symmetry UHF (SEM). Solid line (green): 
EXD (from Ref-H). Thick dashed-dotted line (blue): REM. 
Thin dashed line (violet, marked LLL): the commonly used 
approximate energies £^™ LL (.B) [see Eq. 1131 1. Thin dot- 
ted line (black): E*™ hh (B) (see section II.B). For B < 8 

T, the £ t E S LL (B) and E*™ LL (B) curves coincide; we have 
checked that these curves approach each other also at larger 
values of B, outside the plotted range. Numbers near the bot- 
tom curves denote the value of magic angular momenta \L m , 
see Eq. ©] of the ground state. Corresponding fractional 
filling factors are specified by v = N(N — l)/(2L m ). Param- 
eters used: confinement hu>o — 3.60 meV, dielectric constant 
k = 13.1, effective mass m* = 0.067m e . 



E?£l^{B, L) = E£°{B, L) + V2£: 



EXD 
int. LLL 



(L) — . (13) 



An approximate ground-state energy for the system 
can be found through Eq. (|13|1 by determining the 
angular-momentum value L gs that minimizes this ex- 
pression. In the following, this ground-state energy at a 
given B will be denoted simply by omitting the variable 
L on the left-hand-side of Eq. (|T3J), i.e., Ef^ LL {B) = 

We note that, although few in number, full EXD cal- 
culations for finite B that take into consideration both 
the confinement Huq and the actual complexity of the 
Darwin-Fock spectra (including levels with n > 0) have 
been reported 5,6,8 ^ in the literature for several cases 
with N = 3 and N = 4 electrons. These calculations 
will be of great assistance in evaluating the accuracy of 
the REM method (see section III). 

In the above Eq. (|13f) . we have used exact diagonal- 
ization in the lowest Landau level for evaluating the in- 
terelectron interaction contribution to the total energy. 
In alternative treatments, one may obtain the interelec- 
tron energy contribution through the use of various ap- 
proximate wave functions restricted to the LLL. These 
include the use of the Laughlin wave function and de- 
scendants thereof (e.g., composite fermions), or the ro- 
tating electron wave functions at the limit B — > oo, which 
is reached by setting A = IbV% in the right-hand-side 



of Eq. (defining the displaced orbital). For these 

(B,L), 



cases, we will use the obvious notations E, 
E?*(B,L), and E?™ (B,L). 



Laughlin 
tot, LLL 



III. COMPARISON OF APPROXIMATE 
RESULTS WITH EXACT DIAGONALIZATION 
CALCULATIONS 

A. Ground-state energies in external confinement 

Before proceeding with the presentation of results for 
N > 10, we demonstrate the accuracy of the two-step 
method through comparisons with existing EXD results 
for smaller sizes. In Fig. 1, our calculations for ground- 
state energies as a function of B are compared to EXD 
calculations 5 for N — 4 electrons in an external parabolic 
confinement. The thick dotted line (red) represents the 
broken-symmetry UHF approximation (first step of our 
method), which naturally is a smooth curve lying above 
the EXD one [solid line (green)]. The results obtained 
after restoration of symmetry [dashed-dotted line (blue) ; 
marked as REM] agree very well with the EXD one in 
the whole range 2T<S<15T. 32 We recall here that, 
for the parameters of the QD, the electrons form in the 
intrinsic frame of reference a square about the origin of 
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the dot, i.e., a (0,4) configuration, with the zero indicat- 
ing that no electron is located at the center. According 
to Eq. Lq = 6, and the magic angular momenta are 
given by L m = 6 + 4fc, k = 0, 1, 2, ... 

To further evaluate the accuracy of the two-step 
method, we also display in Fig. 1 [thin dashed line (vio- 
let)] ground-state energies calculated with the commonly 
use d 4 ! 1 ?'? 1 approximation Ef^Yhhi^) ( see section II. B). 
We find that the energies E^^ll(^) t eim to substan- 
tially overestimate the REM (and EXD) energies for 
lower values of B (e.g., by as much as 5.5% at B ~ 4 
T). On the other hand, for higher values of B (> 12 T), 
the energies Ef^^Lhi-B) tend to agree rather well with 
the REM ones. A similar behavior is exhibited also by 
the -Eto?LLL (B) energies [the interaction energies are cal- 
culated within the LLL using the REM wave function; 
dotted line (black)]. We have found that the overesti- 
mation exhibited by the EfJ^Lhi^) energies is due to 
the fact that the actual dimensionless Coulomb coeffi- 
cient £ i nt C LLL(^) [See Eq. fl^jl ] is not independent of the 
magnetic field, but decreases slowly as B decreases when 
the effect of the confinement is considered (see Appendix 
B). A similar agreement between REM and EXD results, 
and a similar inaccurate behavior of the limiting-case ap- 
proximation, was found by us also for N = 3 electrons 
in the range 2 T < B < 16 T shown in Fig. 2 (the EXD 
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FIG. 2: (Color online) Two-step method versus EXD cal- 
culations: Ground-state energies (per particle, referenced to 
hfl) for N = 3 electrons. The electrons are arranged in 
a (0,3) structure in the intrinsic frame of reference. Thick 
dashed line (red): broken-symmetry UHF (SEM). Thinner 
solid line (green): EXD (from Ref. 0). Thick solid line 
(blue): REM. Thin dashed line (violet): the commonly used 
approximate energies -E^lll (B) (see text) . Thin dotted line 
(black): ££uf LL (fl) (see text). For B < 8 T, the E?$> hh {B) 
and lll(-S) curves coincide; we have checked that these 
curves approach each other also at larger values of B, outside 
the plotted range. Numbers denote the value of magic angular 
momenta (L m ) of the ground state. Corresponding fractional 
filling factors are specified by v = N(N — l)/(2L m ). Param- 
eters used: confinement Tiujo = 3.37 meV, dielectric constant 
k — 12 A, effective mass m* = 0.067m e . 



TABLE I: Comparison of yrast-band energies obtained from 
REM and EXD calculations for N = 6 electrons in the low- 
est Landau level, that is in the limit B — *• oo. In this limit 
the external confinement can be neglected and only the in- 
teraction energy contributes to the yrast-band energies. En- 
ergies in units of 6 2 /(kIb)- For the REM results, the (1,5) 
polygonal-ring arrangement was considered. For L < 140, 
see Table IV in Ref. 7(b) and Table III in Ref. @(c). The 
values of the fractional filling may be obtained for each L as 



v = N{N- 


1)/(2L). 






L 


REM 


EXD 


Error (%) 


140 


1.6059 


1.6006 


0.33 


145 


1.5773 


1.5724 


0.31 


150 


1.5502 


1.5455 


0.30 


155 


1.5244 


1.5200 


0.29 


160 


1.4999 


1.4957 


0.28 


165 


1.4765 


1.4726 


0.27 


170 


1.4542 


1.4505 


0.26 


175 


1.4329 


1.4293 


0.25 


180 


1.4125 


1.4091 


0.24 


185 


1.3929 


1.3897 


0.23 


190 


1.3741 


1.3710 


0.23 


195 


1.3561 


1.3531 


0.22 


200 


1.3388 


1.3359 


0.21 


calculation 


was taken from 


Ref. [H). 





In all cases, the total energy of the REM is lower than 
that of the SEM (see, e.g., Figs. 1 and 2). Indeed, a 
theorem discussed in Sec. Ill of Ref. |33l pertaining to 
the energies of projected wave functions, guarantees that 
this lowering of energy applies for all values of N and B. 



B. Yrast rotational band at B — > oo 

As a second accuracy test, we compare in Table I REM 
and EXD results for the interaction energies of the yrast 
band for N = 6 electrons in the lowest Landau level [an 
yrast state is the lowest energy state for a given magic 
angular momentum L m , Eq. @]. The relative error is 
smaller than 0.3%, and it decreases steadily for larger L 
values. 



IV. ILLUSTRATIVE EXAMPLES FROM 
MICROSCOPIC REM CALCULATIONS 

A. Which ring isomer has the lowest ground-state 
energy?: REM versus UHF energies 

For a given number N of electrons, there exis1»2& in 
general more than one polygonal-ring isomers, associ- 
ated with stable and metastable equilibrium configura- 
tions of N electrons inside an external harmonic con- 
finemnet huiQ. Figure 3 displays UHF and REM ground- 
state energies for N = 6 and N — 9 electrons associated 
with the various classical polygonal-ring configurations. 
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FIG. 3: (Color online) Comparison of REM and UHF ground- 
state energies per particle (referenced to hQ) associated with 
different ring isomers for N = 6 and N = 9 electrons as a 
function of the magnetic field B. The curves are labeled with 
the computational method and the isomer (ni, 712). To the left 
of the vertical arrow (at B = 11.5 T), the UHF(1,8) curve is 
energetically favored. To the right of the vertical arrow, the 
UHF(2,7) curve is energetically favored. Parameters used: 
confinement Tiojo = 3.60 meV, dielectric constant n = 13.1, 
effective mass m* = 0.067m e . 



For N — 6, one has two isomers, i.e., a (0,6) configura- 
tion and a (1,5) configuration (with one lectron at the 
center). For N — 9 electrons, there exist three different 
isomers, i.e., (0,9), (1,8), and (2,7). From the bottom 
panel in Fig. 3, we observe that for TV = 9 electrons, the 
lowest REM energies correspond to the classically stable 
isomer, i.e., to the (2,7) configuration with two electrons 
in the inner ring and seven electrons in the outer ring. 
In particular, we note that the (0,9) isomer (which may 
be associated with a single- vortex state) yileds REM en- 
ergies far above the (2,7) one in the whole magnetic-field 
range 5 T < B < 25 T, and in particular for magnetic 
fields immediately above those associated with the MDD 
(the so-called MDD break-up range); the MDD for N = 9 
electrons has an angular momentum Lq = 36 and corre- 
sponds to the first energy oscillation in the figure. 

We have found that the (0, N) isomer is not associated 
with REM ground energies for any magnetic-field range 
in all cases with N > 7. The only instance when the 
(0, N) configuration is associated with a REM ground- 
state energy is the N = 6 case [see Fig. 3, top frame], 
where the REM energy of the (0,6) configuration provides 
the ground-state energy in the range 6.1 T < B < 7.7 T, 
immediately after the break-up of the MDD. 
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FIG. 4: (Color online) Ground-state energies [i.e., for the (2,7) 
configuration] for N = 9 electrons (per particle, referenced to 
HQ) as a function of the magnetic field B. Dashed line (red): 
UHF (SEM). Solid line (blue): REM. Dotted line (black): ap- 
proximate energies Ef^^ Lh (B) (see text). Parameters used: 
confinement fhJo — 3.60 meV, dielectric constant k = 13.1, 
effective mass m* = 0.067m e . 



For comparison, we have also plotted in Fig. 3 the UHF 
energies as a function of the magnetic field. Most notice- 
able is the fact that the REM ground states, compared 
to the UHF ones, may result in a different ordering of 
the isomers. For example, in the range 5 T < B < 6.1 T, 
the UHF indicates, by a small energetic advantage, the 
(0,6) as the ground-state configuration associated with 
the MDD, while the REM specifies the (1,5) arrangement 
as the ground-state configuration. A similar switching of 
the ground-state isomers is also seen between the (1,8) 
and (2,7) configurations in the case of N = 9 electrons in 
the magnetic-field range 5T<£><11.5T. We conclude 
that transitions between the different electron-molecule 
isomers derived from UHF energies alon o 34 i 35 are not re- 
liable. 



B. The case of N = 9 electrons 

In figure 4 we show ground-state energies for the case of 
N = 9 electrons, which have a nontrivial double- ring con- 
figuration (711,712). Here, the most stable configuration 
for classical point charges 25 is (2,7), for which we have 
carried UHF (SEM) and REM (projected) calculations in 
the magnetic field range 5T<£?<25T. We also display 
in Fig. 4 the energies E^™ LL (B) [dotted curve (black)], 
which, as in the N = 4 and N = 3 cases discussed in 
the section III, overestimate the ground-state energies, 
in particular for smaller Bm^ In keeping with the find- 
ings for smaller sizes£( c ) [with (0,N) or (1,N — 1) con- 
figurations], both the UHF and the REM ground-state 
energies of the N = 9 case approach as B — ► 00 the 
classical equilibrium energy of the (2,7) polygonal con- 
figuration [i.e., 16.75 meV; 4.088i?o in the units of Ref. 
M E = (m*w 2 e 4 /2K 2 ) 1/3 ]. 

In analogy with smaller sizes [see, e.g., Figs. 1 and 
2 for N = 4 and N = 3], the REM ground-state 
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TABLE II: Ground-state magic angular momenta and their 
decomposition {ki,k 2 } for N = 9 in the nagnetic- field range 
5 T < B < 25 T. These results correspond to the REM [see 
lower curve in Fig. 4, with the electrons arranged in a (2,7) 
structure]. The parameters used are as in Fig. 4. 
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energies in Fig. 4 exhibit oscillations as a function of 
B. These oscillations reflect the incompressibility of the 
many-body states associated with magic angular mo- 
menta. The magic angular momenta are specified by the 
number of electrons on each ring, and in general they 
are given by L m — N(N — l)/2 + Yjq=i k q n q , where 
the n 9 's are the number of electrons located in the gth 
ring and the fc g 's are non-negative integers; in particu- 
lar, L m = 36 + 2ki + 7k 2 for the N = 9 case in Fig. 4. 
An analysis of the actual values taken by the set of in- 
dices {k\,k 2 } reveals several additional trends that fur- 
ther limit the allowed values of ground-state L m 's. In 
particular, starting with the values {fci = 0, k 2 = 0} at 
B = 5 T (L^ DD = 36), the indices {h,k 2 } reach the 
values {2,24} at B = 25 T (L m = 208). As seen form 
Table II, the outer index k 2 has a short period, while 
the inner index ki exhibits a longer period and increases 
much more slowly than k 2 . This behavior minimizes the 
total kinetic energy of the independently rotating rings 
(having a variable radius, see section V below). 

We also list in Table III the first few pairs of indices 
{ki, k 2 } associated with the £^™ LL curve ( see top dot- 
ted curve in Fig. 4). It can be seen that the magic an- 
gular momenta are different from those associated with 
the REM curve, when the confinement is taken into con- 



TABLE III: Ground-state magic angular momenta and their 
decomposition {fei,fe} for N = 9 electrons associated with 
the -E'tot'.LLL curve [top curve in Fig. 4; see section II. B for an 
explanation of notation; the electrons are arranged in a (2,7) 
structure] . 
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FIG. 5: (Color online) REM radial electron densities for the 
MDD (L m = L = 36) of N = 9 electrons [in the (2,7) ground- 
state configuration] at (a) B — > oo, i.e., in the lowest Landau 
level and (b) at B = 5.5 T. Parameters used in (b): confine- 
ment huio — 3.60 meV, dielectric constant n — 13.1, effective 
mass m* = 0.067m e . Lengths: (a) in units of the magnetic 
length l B ; (b) in units of R = (2e 2 /m* : kwq) 1/3 '. Electron 
densities: (a) in units of 1/1%; (b) in units of 1/Ro- Normal- 
ization: 2n J °° p(r)rdr — N. 



sideration using the full projected wave function in Eq. 
(J2J. The magic angular momenta of the E^^ hL curve 
coincide with the L m 's of the EXD within the LLL, and 
thus are characterized by having Lq + N — 45 (instead 
of Lo + n 2 = 43) as the magic angular momentum imme- 
diately following that of the MDD (i.e., L = 36). The 
L a + N magic angular momentum is associated with the 
(0, N) ring arrangement, which can be interpreted as a 
single "vortex-in-the-center" state. 

Based on EXD calculation^ restricted Jto the lowest 
Landau leveii2i2L2£ (that is, £^ LL or Ef^LL m our 
notation), it has been conjectured that for QDs with 
N < 15, the break-up of the MDD with increasing B 
proceeds through the formation of the above mentioned 
single central vortex state. However, our REM calcu- 
lations show (see also the case of N = 11 electrons in 
section IV. C and the case of N = 17 electrons in section 
IV. D) that taking into account properly the influence of 
the confinement does not support such a scenario. In- 
stead, the break-up of the MDD resembles an edge re- 
construction and it proceeds (for all dots with N > 6) 
through the gradual detachement of the outer ring asso- 
ciated with the classical ground-state polygonal config- 
uration (see Table II for the case of N = 9 electrons). 
The only case we found where the break-up of the MDD 
proceeds via a (0,N) vortex state is the one with N = 6 
electrons (see section IV. A),; and naturally the cases with 
N < 5. 

As another illustration of the subtle, but important, 
differences that exist between wave functions defined ex- 
clusively within the LLL and those specified by the REM 
wave function for finite B in Eq. (0), we display in Fig. 
5 for N = 9 electrons the radial electron densities of the 
MDD at B -> oo and at B = 5.5 T. While the electron 
density of the MDD in the LLL (B — > oo) is constant 
in the central region [up to r » 31b, see Fig. 5(a)], the 
corresponding density at B — 5.5 T displays the charac- 
teristic oscillation corresponding to the (2,7) multi-ring 
structure [see Fig. 5(b)]; the latter behavior is due to the 
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mixing of higher Landau levels. To further illustrate the 
(2,7) crystalline character of the MDD when higher Lan- 
dau levels are considered, we display in Fig. 6 the corre- 
sponding CPDs associated with the REM wave function 
of the MDD at B = 5.5 T and an external confinement 
of hujQ = 3.6 meV. Our conclusions concerning the MDD 
electron densities (and CPDs) are supported by EXD cal- 
culations for N — A electrons^ Note that, while the ring 
structure is well developed in the CPDs shown in Fig. 6, 
the internal (2,7) structure of the rings (see in particu- 
lar the outer ring in the left panel in Fig. 6) is rather 
weak, as expected for the lowest angular momentum Lq 
(MDD). However, the ring structure is easily discernible 
in contrast to the CPDs for the MDD restricted to the 
LLL where structureless CPDs (as well as structureless 
electron densities) are found. 

C. The case of N — 11 electrons 

Figure 7 presents the case for the ground-state energies 
of a QD with N = 11 electrons, which have a nontriv- 
ial double-ring configuration (ni,n 2 ). The most stable^ 
classical configuration is (3, 8), for which we have carried 
UHF (SEM) and REM (projected) calculations in the 
magnetic field range 5T<£?<25T. Figure 7 also dis- 
plays the LLL ground-state energies E^^ LL (B) [dotted 
curve (black)], which, as in previous cases, overestimate 
the ground-state energies for smaller B. The approxi- 
mation £^™ LL (i?), however, can be used to calculate 
ground-state energies for higher values of B. In keep- 
ing with the findings for smaller sizes^- 1 [with (0, N) 
or (1, N — 1) configurations], we found that both the 




FIG. 6: (Color online) Conditional probability distributions 
obtained from REM wave functions of the MDD (Lo = 36) for 
N = 9 electrons at B = 5.5 T [see Fig. 5(b)]. The electrons 
are arranged in a (2,7) structure. The observation point is 
denoted by a solid dot. On the left, the observation point is 
located on the outer shell, and on the right it is located on the 
inner shell. Parameters used: confinement huio — 3.60 meV, 
dielectric constant k = 13.1, effective mass m* = 0.067m e . 
Lengths in units of Ro = (2e 2 /(/tm*Wo)) 1 ■ CPDs (vertical 
axes) in arbitrary units. 
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FIG. 7: (Color online) Ground-state energies for N = 11 
electrons (per particle, referenced to hQ.) as a function of 
the magnetic field B. Dashed line (red): UHF (SEM). Solid 
line (blue): REM. Dotted line (black): Approximate ener- 
gies EtetjuLh(B) (see text). Parameters used: confinement 
Hljo = 3.60 meV, dielectric constant n = 13.1, effective mass 
m* = 0.067m e . The inset shows a magnification of the REM 
curve in the range 5 T < B < 12 T. 

UHF and the REM ground-state energies approach, as 
B — > oo, the classical equilibrium energy of the (3,8) 
polygonal co nfig uration [i.e., 19.94 meV; 4.865-Eo in the 
units of Ref.Hl E = (m*^e 4 /2k 2 ) 1/3 ]- 

In analogy with smaller sizes [see, e.g., Figs. 1, 2, and 
4 for TV = 4, 3, and 9, respectively], the REM ground- 
state energies in Fig. 7 exhibit oscillations as a function 
of B (see in particular the inset) . As discussed in section 
IV. B, these oscillations are associated with magic angular 
momenta, specified by the number of electrons on each 
ring. For N — 11 they are given by Eq. ©, i.e., L m — 
55 + 3/ji + 8fc2, with the kq's being nonnegative integers. 
As was the case with N = 9 electrons, an analysis of the 
actual values taken by the set of indices {fci, A^} reveals 
several additional trends that further limit the allowed 
values of ground-state L m 's. In particular, starting with 



TABLE IV: Ground-state magic angular momenta and their 
decomposition {fci, £2} for N — 11 in the nagnetic-field range 
5 T < B < 25 T. The results correspond to the REM (see 
lower curve in Fig. 7). The parameters used are as in Fig. 7. 
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FIG. 8: (Color online) REM conditional probability distribu- 
tions for N = 11 electrons at B = 10 T (L = 106). The 
electrons are arranged in a (3,8) structure. The observa- 
tion point (solid dot) is placed on (left) the outer ring at 
ro = 1.480i?o, and (right) on the inner ring at ro = 0.557i?o- 
Parameters used: confinement hu>o — 3.60 meV, dielectric 
constant k = 13.1, effective mass m* = 0.067m,;. Lengths 
in units of Ro = (2e 2 /m *kujI) 1/3 . CPDs (vertical axes) in 
arbitrary units. 



the values {0,0} at B = 5 T (L = 55), the indices 
{fci,fc 2 } reach the values {3, 24} at B = 25 T (L m = 256). 
As seen from Table IV, the outer index /c2 changes faster 
than the inner index k\. This behavior minimizes the 
total kinetic energy of the independently rotating rings; 
indeed, the kinetic energy of the inner ring (as a function 
of fci) rises faster than that of the outer ring (as a function 
of fo) due to smaller moment of inertia (smaller radius) 
of the inner ring [see Eq. I|14l) ]. 

In addition to the overestimation of the ground-state 
energy values, particularly for smaller magnetic fields 
(see Fig. 7 and our above discussion), the shortcomings 
of the LLL approximation pertaining to the ground-state 
ring configurations [see section II. B, Eq. I|13fl ]. as dis- 
cussed by us above for N = 9, persist also for N = 11. 
In particular, we find that according to the LLL ap- 
proximation the ground-state angular-momentum imme- 
diately after the MDD (Lq = 55) is L m — 66, i.e., the 
one associated with the (0, N) vortex-in-the-center con- 
figuration. This result, erroneously stated in Ref. as 
the ground state, disagrees with the correct result that 
includes the effect of the confinement - listed in Table IV, 
where the ground-state angular momentum immediately 
following the MDD is L m = 63. This angular momen- 
tum corresponds to the classicaly most stable (3,8) ring 
configuration - that is a configuration with no vortex at 
all. 

Figure 8 displays the REM conditional probability dis- 
tributions for the ground state of TV = 11 electrons at 
B = 10 T (L m = 106). The (3,8) ring configuration 
is clearly visible. We note that when the observation 
point is placed on the outer ring (left panel), the CPD 
reveals the crystalline structure of this ring only; the in- 



ner ring appears to have a uniform density. To reveal the 
crystalline structure of the inner ring, the observation 
point must be placed on this ring; then the outer ring 
appears to be uniform in density. This behavior suggests 
that the two rings rotate independently of each other, a 
property that we will explore in section V to derive an 
approximate expression for the yrast rotational spectra 
associated with an arbitrary number of electrons. 



D. The case of N = 17 electrons 

Figure 9 presents (for 5 T < B < 15 T) REM and 
UHF ground-state energies for TV = 17 electrons, which 
have a (1,6,10) three-ring configuration as the most stable 
classical arrangement^ 

In analogy with smaller sizes [see, e.g., previous figures 
for N < 12] the REM ground-state energies in Fig. 9 ex- 
hibit oscillations as a function of B, and each oscillation 
is associated with a given particular (magic) value of the 
angular momentum. Earlier in this section we discussed 
the physical origins of the magic angular momenta. As 
before, the magic angular momenta are specified by the 
number of electrons on each ring [0], with L = 136 
and L m = 136 + 6fc 2 + 10fc 3 for N = 17; fc g 's being 
non-negative integers (the central electron does not con- 
tribute to the total angular momentum). Analysis of 
the particular values taken by the set of indices {/c2, k^} 
reveals similar trends to those found for the cases with 
N = 9 and TV = 11 electrons. In particular, starting 
with the values {0, 0} at B = 5 T (Lq = 136), the indices 
{&2, fe} reach the values {&2 = 5, k 3 = 18} at B = 15 T 
(L m = 346). As seen from Table V, the outer index k 3 
changes faster, than the inner index k^- This behavior 
minimizes the total kinetic energy of the independently 
rotating rings, as was already discussed for iV = 9 and 
N — 11 electrons. 
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FIG. 9: (Color online) Ground-state energies (per particle, 
referenced to hQ) for N — 17 electrons as a function of the 
magnetic field B. The electrons are arranged in a (1,6,10) 
structure. Dashed line (red): UHF. Solid line (blue): REM. 
Parameters used: confinement Tiuoo — 3.6 meV, dielectric con- 
stant « = 13.1, effective mass m* = 0.067m e . 
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TABLE V: Ground-state magic angular momenta and their 
decomposition {^2,^3} for N = 17 electrons in the nagnetic- 
field range 5 T < B < 15 T. The rersults correspond to the 
REM (see lower curve in Fig. 9) . The parameters used are as 



in Fig. 9. 
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FIG. 11: (Color online) CPDs for N = 12 electrons and with 
angular momentum L = 132 (u — 1/2) calculated with EXD 
in the lowest Landau level. The electrons are arranged in a 
(3,9) structure. The observation point (solid dot) is placed on 
the outer ring at ro = 5.221b (left frame), and on the inner 
ring at ro = 1.871b (right frame). Lengths in units of Is- 
CPDs (vertical axes) in arbitrary units. 



We have also calculated the ground-state energies for 
N = 17 electrons in the LLL approximation, i.e., we cal- 
culated the quantity E^^ hh (B) (not shown in Fig. 9). 

We find once more that E^^ LL (B) overestimates the 
ground-state energies in the magnetic-field range cov- 
ering the MDD and the range immediately above the 
MDD. For N = 17, however, the shortcoming of the LLL 
approximation is not reflected in the determination of 
the ground-state ring configurations. We find that for 
N = 17 the LLL approximation yields a (1,6,10) ring con- 
figuration (with L rn — 146) for the ground state imme- 
diately following the MDD, in agreement with the REM 
configurations listed in Table V. 




FIG. 10: (Color online) Grond-state conditional probability 
distributions, CPDs, obtained from REM wave functions for 
the ground state of TV = 17 electrons at B = 10 T (L = 
228). The electrons are arranged in a (1,6,10) structure. The 
observation point (solid dot) is placed on the outer ring at 
ro = 1.858-Ro (left frame), and on the inner ring at ro = 
0.969i?o (right frame). The rest of the parameters are the 
same as in Fig. 9. Lengths in units of Ro = (2e 2 /(/tm*Uo)) 1 ^ 3 - 
CPDs (vertical axes) in arbitrary units. 



V. REM YRAST BAND EXCITATION 
SPECTRA AND DERIVATION OF ANALYTIC 
APPROXIMATE FORMULA 

In Fig. 10, we display the CPD for the REM wave 
function of N = 17 electrons. This case has a nontriv- 
ial three-ring structure (1,6,10), 2 - which is sufficiently 
complex to allow generalizations for larger numbers of 
particles. The remarkable combined character (partly 
crystalline and partly fluid leading to a non-classical ro- 
tational inertia, see section VI) of the REM is illustrated 
in the CPDs of Fig. 10. Indeed, as the two CPDs [re- 
flecting the choice of taking the observation point [ro 
in Eq. JSj] on the outer (left frame) or the inner ring 
(right frame)] reveal, the polygonal electron rings rotate 
independently of each other. Thus, e.g., to an observer 
located on the inner ring, the outer ring will appear as 
having a uniform density, and vice versa. The wave func- 
tions obtained from exact diagonalization exhibit also 
the property of independently rotating rings [see e.g., the 
N = 12 and L = 132 {v = 1/2) case in Fig. 11], which is 
a testimony to the ability of the REM wave function to 
capture the essential physics^ of a finite number of elec- 
trons in high B. In particular, the conditional probabil- 
ity distribution obtained for exact diagonalization wave 
functions in Fig. 11 exhibits the characteristics expected 
from the CPD evaluated using REM wave functions for 
the (3,9) configuration and with an angular-momentum 
decomposition into shell contributions [see Eqs. J2J an d 
0] La = 3 + 3fci and L 2 = 63 + 9k 2 {L x + L 2 = L m ; 
for L m = 132 the angular-momentum decomposition is 
Li = 6 and L 2 = 126). 

In addition to the conditional probabilities, the 
solid/fluid character of the REM is revealed in its ex- 
cited rotational spectrum for a given B. From our micro- 
scopic calculations based on the wave function in Eq. J2J) , 
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we have derived (see below) an approximate (denoted as 
"app"), but analytic and parameter- free, expression [see 
Eq. 119J1 below] which reflects directly the nonrigid (non- 
classical) character of the REM for arbitrary size. This 
expression allows calculation of the energies of REMs for 
arbitrary N, given the corresponding equilibrium config- 
uration of confined classical point charges. 

We focus on the description of the yrast band at a 
given B. Motivated by the aforementioned nonrigid 
character of the rotating electron molecule, we consider 
the following kinetic-energy term corresponding to a 



(m, . 



.,n r ) configuration (with . n q — N): 



h 2 L 2 /(2J q (a q )) - hco c L/2, (14) 



where L q is the partial angular momentum associated 
with the qth ring about the center of the dot and the total 
angular momentum is L = J2 q =i Eq- J q {a q )) = n q m*a q 
is the rotational moment of inertia of each individual 
ring, i.e., the moment of inertia of n q classical point 
charges on the gth polygonal ring of radius a q . To ob- 
tain the total energy, Ef EM , we include also the term 
E *p P ( N ) = Eq=i J q {a q )n 2 /2 due to the effective har- 
monic confinement f2 [see discussion of Eq. QJ], as well 
as the interaction energy E£ , 

r Q 2 r— 1 r 



q=l 



4 kq 



(15) 



q=l s>q 



The first term is the intra-ring Coulomb-repulsion energy 
of n q point-like electrons on a given ring, with a structure 
factor 



j=2 



(sin[(j - l)7r/n g ])- 



(16) 



The second term is the inter-ring Coulomb-repulsion en- 
ergy between rings of uniform charge distribution cor- 
responding to the specified numbers of electrons on the 
polygonal rings. The expression fo Vc is 

V c {a q ,a s ) = n q n s 2 F 1 [3/4,l/4;l;Aa 2 q a 2 s (a 2 q + a 2 s y 2 } 



xe*(al+al)-^/ K , 



(17) 



where 2-F1 is the hypergeometric function. 

For large L (and/or B), the radii of the rings of the 
rotating molecule can be found by neglecting the inter- 
action term in the total approximate energy, thus mini- 
mizing only E*™(N) + E^ p (N). One finds 



= A 



\jL q /n q ; 



(18) 



i.e., the ring radii depend on the partial angular momen- 
tum L q , reflecting the lack of radial rigidity. Substitution 
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FIG. 12: (Color online) Left: Yrast spectrum for N = 17 
electrons at a high magnetic field B = 100 T. Approxi- 
mate analytic expression [Eq. 1191 . dashed line (violet)] com- 
pared with microscopic REM calculations [Eq. (jj), solid line 
(green)]. Right: The corresponding classical (rigid rotor) en- 
ergy E T ^ S for N — 17 electrons (see text). The microscopic 
REM energies are referenced relative to the zero-point energy, 
17RS1. Energies were calculated for magic angular momenta 
L = Li + L 2 + L 3 with Li = 0, Li = 21 + 6fa> and fc 2 = 30, 
and L3 = 115+ IO/C3. The parameters are the same as in Fig. 
9. Note the much larger energy scale for the classical case 
(right frame), leading to a non-rigidity index for the REM of 
a ~ 0.99 (see text). 



into the above expressions for i?app > -^app > an d -^app yields 
for the total approximate energy the final expression: 



£ a R pp M L {N) = niSl-u c /2)L- 



q=l -^q q= \ s>q 

where the constants 

C v , q = Q.2bn q /2 S q e 2 /{n\). 



(20) 



For simpler (0, N) and (1, TV— 1) ring configurations, Eq. 
(|19fl reduces to the expressions reported earlier. 7 ( c )di 



VI. A NON-RIGID CRYSTALLINE PHASE: 
NON-CLASSICAL ROTATIONAL INERTIA OF 
ELECTRONS IN QUANTUM DOTS 

In Fig. 12 (left frame), and for a sufficiently high mag- 
netic field (e.g., B — 100 T such that the Hilbert space 
of the system reduces to the lowest Landau level), we 
compare the approximate analytic energies E^^ with 
the microscopic energies Ef EM calculated from Eq. © 
using the same parameters as in Fig. 9 for N — 17 elec- 
trons. The two calculations agree well, with a typical 
difference of less than 0.5% between them. More impor- 
tant is the marked difference between these results and 
the total energies of the classical (rigid rotor) molecule 
(E^ lg ), plotted in the right frame of Fig. 12; the latter 
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are given by 

N 

E L ig = h 2 L 2 /(2J lig ) + 0.5j2 m * UJ l\ z ^\ 2 + 

»=i 

N N 

Y J Y, e Vm-z j \), (2i) 

i—l j>i 

with the rigid moment of inertia being 42 

N 

J lig = J2 m *\ Z i\ 2 - ( 22 ) 
i=l 

The disagreement between the REM and the classical 
energies is twofold: (i) The L dependence is different, 
and (ii) The REM energies are three orders of magnitude 
smaller than the classical ones. That is, the energy cost 
for the rotation of the REM is drastically smaller than 
for the classical rotation, thus exhibiting non-classical 
rotational behavior. In analogy with Ref. 18, we define a 
"non-rigidity" index 

a = (E L iB - Ef EM )/E L is . (23) 

For the case displayed in Fig. 12, we find that this in- 
dex varies (for 1116 < L < 3716) from a = 0.978 to 
a = 0.998, indicating that the rotating electron molecule, 
while possessing crystalline correlations is (rotationally) 
of a high non-rigid nature. We remark that our definition 
of a in Eq. I|23|) was motivated by a similar form of an in- 
dex of supefluid fraction introduced in Ref. ll8t we do not 
mean to imply the presence of a superfluid component 
for electrons in quantum dots. 

In the context of the appearance of supersolid behav- 
ior of 4 He under appropriate conditions, formation of a 
supersolid fraction is often discussed in conjunction with 
the presence of (i) real defects and (ii) real vacancies iILii 
Our REM wave function [Eq. [5] belongs to a third possi- 
bility, namely to virtual defects and vacancies, with the 
number of particles equal to the number of lattice sites 
(in the context of 4 He, the possibility of a supersolid with 
equal number of particles and lattice sites is mentioned in 
Ref. l20f) . Indeed, the azimuthal shift of the electrons by 
(T1i72) •••j7r) [ see Eq. may be viewed as generating 
virtual defects and vacancies with respect to the original 
electron positions at (71 =0,72 = 0, ...,7 r = 0) on the 
polygonal rings. 

A recent publication^! has explored the quantal nature 
of the 2D electron molecules in the lowest Landau level 
(B — > 00) using a modification of the second-quantized 
LLL form of the REM wave functions^ In particular, the 
modification consisted of a multiplication of the param- 
eter free REM wave function by variationally adjustable 
Jastrow-factor vortices. Without consideration of the ro- 
tational properties of the modified wave function, the 
inherently quantal nature of the molecule was attributed 
exclusively to the Jastrow factor. However, as shown 
above, the original REM wave function [Eq. J2J] already 



exhibits a characteristic non-classical rotational inertia 
(NCRI). Consequently, the additional variational free- 
dom introduced by the Jastrow prefactor may well lead 
energetically to a slight numerical improvement, but it 
does not underlie the essential quantal physics of the sys- 
tem. Indeed, as discussed previously and illustrated in 
detail above, the important step is the projection of the 
static electron molecule onto a state with good total an- 
gular momentum [see Eqs. and (J2J]. 



VII. SUMMARY 

The focus of this study pertains to the development 
of methods that permit investigations of the energetic, 
structural, and excitation properties of quantum dots 
in strong magnetic fields with an (essentially) arbitrary 
number of electrons. Towards this aim, we utilized sev- 
eral computational methods, and have assessed their ad- 
equacy. The methods that we have used are: (1) Exact 
diagonalization which is limited to a rather small number 
of particles; (2) The "two-step" successive-hierarchical- 
approximations method (see section II. A), in which a 
UHF step leading to broken-symmetry solutions (static 
electron molecule) is followed by restoration (via projec- 
tion techniques) of circularly symmetric states with good 
angular momenta (rotating electron molecule; REM); (3) 
An approximation method based on diagonalization of 
the electron-electron interaction term restricted to the 
lowest Landau level (LLL). In this method, the total 
energy includes, in addition to the LLL diagonalization 
term, a contribution from the harmonic confinement that 
is linear in the total angular momentum; (4) An analytic 
expression [see section V, Eq. Ijl9|l ] whose derivation is 
based on the REM. 

We performed comparative calculations for quantum 
dots with an increasing number of parabolically confined 
electrons (N = 3, 4, 6, 9, 11, and 17). The ground- 
state arrangements of the electrons become structurally 
more complex as the number of electrons in the dot in- 
creases. Using the notation (ni, 72,2, 72.3, ...) for the number 
of electrons located on concentric polygonal rings (see 
section II. A), the ground-state arrangements are: (0,3) 
for N = 3, (0,4) for N = 4, (1,5) for N = 6, (2,7) for 
N = 9, (3,8) for N = 11, and (1,6,10) for N = 17. 

Analysis of the results of our calculations revealed that, 
for all sizes studied by us, the two-step REM method pro- 
vides a highly accurate description of electrons parabol- 
ically confined in quantum dots for a whole range of ap- 
plied magnetic fields, starting from the neighborhood of 
the so-called maximum density droplet and extending to 
the B — > 00 limit. In contrast, the LLL-diagonalization 
approximation was found to be rather inaccurate for 
weaker magnetic fields, where it grossly overestimates the 
total energies of the electrons; the accuracy of this latter 
method improves at higher field strengths. 

The ground-state energy of the electrons in a QD os- 
cillates as a function of the applied magnetic field, and 
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the allowed values of the angular momenta are limited to 
a set of magic angular momentum values, L m , which are 
a natural consequence of the geometrical arrangement of 
the electrons in the rotating electron molecule. Accord- 
ingly, the electrons are localized on concentric polygo- 
nal rings which rotate independently of each other (as 
observed from the conditional probability distributions, 
see section IV). Underlying the aforementioned oscilla- 
tory behavior is the incompressibility of the many-body 
states associated with the magic angular momenta. The 
general expression for L m is given in Eq. @, for a given 
number N and occupancy of the polygonal rings {n q }- 
For the ground-state L m 's, the values of the non- negative 
integers k q in Eq. © are taken such as to minimize the 
total kinetic energy of the electrons. Since the moment 
of inertia of an outer ring is larger than that of an in- 
ner ring of smaller radius, the rotational energy of the 
outer ring will increase more slowly with increasing an- 
gular momentum. Therefore, the k q index in Eq. of an 
outer ring will very up to relatively large values while the 
values corresponding to inner rings remain small (see sec- 
tion IV). As a consequence, we find (see section IV. B to 
IV. D) through REM calculations with proper treatment 
of the confining potential that for N > 6, with increas- 
ing strength of the magnetic field, the maximum density 
droplet converts into states with no central vortex, in 
contrast to earlier conclusionsiSi^LS drawn on the basis 
of approximate calculations restricted to the lowest Lan- 
dau level. Instead we find that the break-up of the MDD 
with increasing B proceeds through the gradual detach- 
ment of the outer ring associated with the corresponding 
classical polygonal configuration. 

In addition to the ground-state geometric arrange- 
ments, we have studied for certain sizes higher-energy 
structural isomers (see, e.g., the cases of N = 6 and 
N = 9 confined electrons in Fig. 3). We find that for 
all cases with N > 7 multi-ring confined-electron struc- 
tures (rii, ri2, n r ), with nx,U2, ...,n r ^ and r > 2, 
are energetically favored. For N — 6, a (1,5) structure is 
favored except for a small £?-range (e.g., 6.1 T < B < 7.7 
T for the parameters in Fig. 3), where the (0,6) single- 
ring structure is favored. For N < 5 the (0, N) single-ring 
structure is favored for all B values. 

In the REM calculations, we have utilized an analytic 
many-body wave function [Eq. J5J] which allowed us to 
carry out computations for a sufficiently large number of 
electrons (N = 17 electrons having a nontrivial three- 
ring polygonal structure), leading to the derivation and 
validation of an analytic expression Eq. Ijl9|l for the total 
energy of rotating electron crystallites of arbitrary N. 

The non-rigidity implied by the aforementioned inde- 
pendent rotations of the individual concentric polygonal 
rings motivated us to quantify (see section VI) the 
degree of non-rigidity of the rotating electron molecules 
at high B, in analogy with the concept of non-classical 
rotational inertia used in the analysisi&SS of supersolid 
4 He. These findings for finite dots suggest a strong 
quantal nature for the extended Wigner crystal in the 



lowest Landau level, designating it as a useful paradigm 
for exotic quantum solids. 
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APPENDIX A: PROOF THAT u(z, Z) [EQ. ©] 
LIES IN THE LLL WHEN A = IbV% 

Using the identity — i{xY - yX) = (zZ* — z*Z)/2, one 
finds 



u(z, Z; A = IbV%) 

zz* -ZZ* +2zZ* 
4/ ? 

e b 



'2-kIb V2ttI b 

OO 

£ci(z*),M*), 



1 zZ* 



11 V 2/| 



1=0 



where z = x + iy, Z = X + iY, and 



with 



Ci(Z*) = 



= ±f. 



\ibV2 



' -ZZ* 

771 



IB V.B 



V2 



e *'b 



(Al) 



(A2) 



(A3) 



being the Darwin-Fock single-particle wave functions 
with zero nodes forming the LLL. 



APPENDIX B: COULOMB MATRIX ELEMENTS 
BETWEEN DISPLACED GAUSSIANS [EQ. ©] 

We give here the analytic expression for the Coulomb 
matrix elements, 

Vijki= [ [ dr 1 dr 2 u*(ri)u*(r 2 )- r ^ r u fe (r 1 )w ; (r 2 ), 

J J k\ti - r 2 | 

(Bl) 

between displaced Gaussians [see Eq. QJ] centered at four 
arbitrary points Zi, Zj, Z^, and Z\. 
One has 



(B2) 



with 



Z t Z* + Z,Z* + Z k Zl + ZiZT 
<? = -— ' 3 2A2 L± + C?; + (TT , (B3 ) 
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and 



(B4) 



The magnetic-field dependence is expressed through 
the parameter 



where 



Zl- ~\~ Z-t Zl? Zn 

c = k J" - +p- 



2A 



2A 



2A 



2A 



(B5) 



(B6) 



\2 



(B9) 



Zi + Zj + Z t - Zj 



2A 



2A 



2A 



Z i+ Z i , a Z i~ Z i 



2A 



The length parameters A and Ib (magnetic length) are 
(B7) defined in the text following Eq. Q. Note that j3 = 
for B = and /3 = 1 for _B — > oo. The latter offers an 
alternative way for calculating REM energies and wave 
functions in the lowest Landau level without using the 
(B8) analytic REM wave functions presented in Ref. \& 
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